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Electrified  vehicles  (EV)  and  renewable  power  sources  are  two  important  technologies  for  sustainable 
ground  transportation.  If  left  unmitigated,  the  additional  electric  load  could  over-burden  the  electric 
grid.  Meanwhile,  a  challenge  for  integrating  renewable  power  sources  into  the  grid  lies  in  the  fact  their 
intermittency  requires  more  regulation  services  which  makes  them  expensive  to  deploy.  Fortunately, 
EVs  are  controllable  loads  and  the  charging  process  can  be  interrupted.  This  flexibility  makes  it  possible 
to  manipulate  EV  charging  to  reduce  the  additional  electric  load  and  accommodate  the  intermittency  of 
renewable  power  sources.  To  illustrate  this  potential,  a  two-level  optimal  charging  algorithm  is  designed, 
which  achieves  both  load  shifting  and  frequency  regulation.  Load  shifting  can  be  realized  through  coor¬ 
dination  of  power  generation  and  vehicle  charging  while  reducing  power  generation  cost  and  carbon 
dioxide  emissions.  To  ensure  practicality,  a  decentralized  charging  algorithm  for  load  shifting  is  formu¬ 
lated  by  emulating  the  charging  pattern  identified  through  linear  programming  optimization  solutions. 
The  frequency  regulation  is  also  designed  based  on  frequency  droop  that  can  be  implemented  in  a  decen¬ 
tralized  way.  The  two  control  objectives  can  be  integrated  because  they  are  functionally  separated  by  time 
scale.  Simulation  results  are  presented  to  demonstrate  the  performance  of  the  proposed  decentralized 
algorithm. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Electrified  vehicles  and  renewable  power  sources  are  actively 
pursued  as  clean  energy  technologies.  However,  the  former  can 
increase  the  electric  grid  load  during  charging,  which  could  stress 
generation  and  transmission  systems.  The  latter  are  intermittent  in 
nature  and  can  bring  challenging  regulation  problems  when  they 
are  integrated  with  the  grid. 

In  the  past,  increased  electric  demand  has  mostly  been  met  by 
building  extra  centralized  power  plants.  The  electric  infrastruc¬ 
ture  was  designed  to  meet  the  peak  demands  which  only  occur 
a  few  hundred  hours  each  year  in  the  US  [1  ].  The  peak  power  costs 
more  because  it  uses  more  expensive  fuels  and  have  shorter  run¬ 
ning  hours  than  the  base  powers  do.  If  the  peak  demand  can  be 
redistributed  to  less  congested  hours,  such  as  during  early  morn¬ 
ing,  the  additional  plants  for  the  peak  power  might  not  need  to  be 
built.  This  is  the  so-called  load  shifting  strategy.  Electrified  vehi¬ 
cles  are  controllable  loads  and  can  even  be  power  sources  (i.e.,  they 
can  return  energy  back  to  the  grid)  under  extraneous  situations  [2]. 
Most  EVs  are  plugged-in  for  long  hours  during  the  night,  i.e.,  exactly 
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during  the  valley  hours  of  traditional  electric  grid.  Therefore,  we  can 
manipulate  the  charging  pattern  of  EVs  to  achieve  load  shifting  or 
valley  filling. 

Regulation  is  a  short  time  scale  ancillary  service  usually  con¬ 
ducted  by  fast-reacting  power  sources.  The  power  sources  for 
regulation  respond  to  a  system-operator’s  request  to  track  the 
minute-by-minute  fluctuations  in  system  load  and  to  respond 
to  unintended  fluctuations  in  generator  output  [3].  They  are 
expensive  to  operate  due  to  the  stand-by  characteristics  of  the 
service  and  their  responsiveness.  When  wind  and  solar  power 
sources  are  introduced,  they  require  a  higher  level  of  regulation 
services,  resulting  in  increased  regulation  cost.  One  of  the  benefits 
of  EVs  is  their  response  time.  Due  to  this  characteristic,  EVs  can 
replace  some  of  the  regulation  service  units. 

Different  approaches  to  realizing  these  two  potential  benefits 
(load  shifting  and  regulation)  have  been  presented  in  the  litera¬ 
ture.  The  cases  in  which  vehicles  provide  power  to  the  grid,  V2G 
(vehicle-to-grid),  have  been  studied  [4,5].  A  number  of  studies  [6-9] 
focused  on  the  regulation  capability  of  EV  batteries  and  have  shown 
that  V2G  may  be  beneficial  to  the  grid  operations,  sometimes  at  the 
expense  of  lower  battery  charging  completion.  Furthermore,  these 
approaches  are  based  on  a  short  time  horizon  and  have  not  explored 
long-term  behaviors  such  as  valley-filling  or  load  shifting.  Even 
though  the  benefit  of  EV  load  shifting  has  been  studied  in  many 
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Nomenclature 

Acronyms 

EV 

electrified  vehicle 

V2G 

vehicle  to  grid 

AC 

alternating  current 

SOC 

state  of  charge 

LP 

linear  programming 

Symbols 

P 

power 

r2,r3 

time  constants  of  power  plants 

ku  kP 

integral  and  proportional  gains  of  a  controller 

CO 

frequency 

H 

normalized  inertial  constant 

D 

normalized  damping  constant 

Pev 

charging  power  of  EV 

PeV max 

the  maximum  charging  power  of  EV 

P nonEV 

non-EV  power  demand 

Pev,sft 

EV  charging  power  that  contributes  demand  shifting 

Pev, reg 

EV  charging  power  that  contributes  frequency  reg¬ 
ulation 

Pev, reg 

EV  charging  power  that  contributes  frequency  reg¬ 
ulation 

NEv 

the  number  of  EVs  connected  to  the  grid 

C(n) 

battery  capacity  of  EV  #n 

papers  [1,10-12],  few  papers  have  presented  how  to  realize  EV 
load  shifting.  Recently,  some  papers  [13,14]  have  begun  to  design 
non-centralized  strategies  for  valley-filling.  Ma  et  al.  [13]  and  Call¬ 
away  and  Hiskens  [14]  showed  that  a  demand  dependent  pricing 
scheme  drives  a  unique  Nash  equilibrium  that  results  in  valley¬ 
filling.  This  algorithm  works  in  a  decentralized  way  but  requires 
that  all  vehicles  have  access  to  information  of  other  vehicles  and  the 
grid  power  generation,  so  the  developed  algorithm  requires  a  lot  of 
communications.  Furthermore,  it  does  not  provide  a  closed  form  of 
the  charging  algorithm.  In  theory,  a  centralized  controller  can  col¬ 
lect  full  information  of  all  EVs  and  all  power  plants,  utilizing  future 
power  demand  and  control  all  vehicles/plants  simultaneously  for 
optimal  performance.  Such  an  approach,  however,  requires  exten¬ 
sive  bi-directional  communication  and  heavy  computations  and 
thus  is  not  desirable.  A  decentralized  charging  controller  requires 
smaller  amount  of  information  to  be  collected  and  exchanged. 
Most  of  the  existing  algorithms  of  load  shifting  are  centralized  and 
thus  require  more  infrastructure  for  implementation.  Furthermore, 
most  of  them  are  based  on  thermal  load  [  1 5,1 6]  or  air-conditioning 
load  [17]  control  not  on  EV  load  control.  Besides,  most  papers 
[6-9,13,14]  on  the  V2G  framework  presented  only  one  control 
objective,  regulation  or  load  shifting,  but  not  both. 

In  our  design,  the  decentralized  and  distributed  charging  con¬ 
trollers  receive  a  simple  command  from  the  centralized  grid 
controller  and  they  only  have  access  to  the  information  from  the 
local  vehicle  (e.g.,  battery  state  of  charge)  and  the  local  grid  (e.g., 
AC  frequency).  Designing  such  decentralized  charging  control  algo¬ 
rithms  that  achieve  near-optimal  performance  (compared  to  a 
centralized  controller)  is  the  goal  and  the  main  contribution  of  this 
paper.  In  this  paper,  we  present  a  decentralized  charging  control 
algorithm  that  minimizes  electricity  generation  cost  and  carbon 
dioxide  emissions.  The  algorithm  also  regulates  charging  power  to 
mitigate  uncertainties  due  to  renewable  power  sources.  The  con¬ 
troller  is  developed  with  two  separate  time  scales:  a  load  shifting 
algorithm  on  a  long  time  scale  and  a  frequency  regulation  algorithm 
on  a  short  time  scale.  Effects  of  the  two  algorithms  are  sepa¬ 
rated  due  to  time  scale  and  the  algorithms  can  be  easily  combined 


without  compromising  their  performance.  The  remainder  of  the 
paper  is  organized  as  follows:  the  problem  is  defined  in  Section  2; 
the  design  of  the  load  shifting  algorithm  is  explained  in  Section  3; 
the  frequency  regulation  algorithm  is  presented  in  Section  4;  the 
integration  and  overall  structure  and  performance  of  the  controller 
are  presented  in  Sections  5  and  6;  and  finally  a  conclusion  is  given 
in  Section  7. 

2.  Problem  statement 

2.1.  Grid  and  demand 

In  this  paper,  the  area  supported  by  DTE  Energy  in  Michigan  is 
used  to  determine  the  size  the  power  generation  and  load  [18].  The 
electric  power  demand  in  this  area  varies  between  5000  MW  and 
8000  MW,  with  the  summer  load  higher  than  the  winter  load  due  to 
air-conditioning.  Peak  electric  demand  occurs  around  2  pm  in  the 
summer  and  7  pm  in  the  winter.  The  lowest  demand  occurs  around 
4  am.  The  hourly  power  demand  profiles  for  typical  summer  and 
winter  days  are  shown  in  Fig.  1,  where  the  load  below  5500  MW 
is  defined  as  the  base  load,  the  load  between  5500  and  7200  MW 
is  defined  as  intermediate  load,  and  above  7200  MW  is  defined  as 
peak  load.  The  summer  profile  will  be  used  for  this  paper.  The  same 
design  process  of  course  can  be  used  for  the  winter  load. 

2.2.  Generation  units  and  frequency  dynamics 

All  electric  grid  systems  have  multiple  types  of  generation 
sources.  The  first  type  is  base  load  sources  that  have  the  cheapest 
operating  cost.  The  second  type  is  load  following  units  which  are 
used  to  make  up  the  difference  between  the  demand  and  the  base 
load  power.  The  third  type  of  power  sources  are  regulation  units 
and  are  used  for  frequency  regulation.  The  characteristics  of  the 
three  types  of  power  sources  are  summarized  in  Table  1.  The  base 
load  power  plants  are  assumed  to  follow  a  planned  profile  and  be 
the  same  as  the  forecasted  power  demand.  The  load  following  units 
and  regulation  units  are  modeled  as  first  order  systems  controlled 
by  an  integral  controller  and  a  proportional  controller,  respectively 
[19],  as  shown  in  Eqs.  (1)  and  (2) 

^1  —  ^planned)  P 2  —  ^  _|_  ~^/oIlow»  ^3  —  ^  ^“^freg>  (1) 

where  Pi ,  P2  and  P3  are  the  power  output  from  the  three  types 
of  power  source.  r2  and  r3  are  the  time  constants  for  the  second 
and  third  type  of  power  plants.  The  power  generation  strategy  to 
regulate  the  frequency  is  simplified  as  follows: 

Pfoilow  —  ~ki  J  A co  dt ,  Preg  =  —kpAco,  (2) 

where  1q  and  kp  are  control  gains  designed  to  stabilized  the  system 
and  A  co  is  a  deviation  from  the  desired  frequency  value,  60  Hz  in 
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Fig.  2.  Power  demand  profile  when  the  charging  of  the  2  M  EVs  are  unmitigated 
(charging  at  the  maximum  power  level  immediately  after  plug  in). 

US.  The  frequency  of  the  grid  is  modeled  as  first-order  rotational 
dynamics,  as  follows: 

Hco  +  Deo  =  Pm  -  Pe ,  (3) 

where  H  is  the  inertial  constant  in  normalized  (per)  unit,  D  is  a 
damping  coefficient  in  per  unit,  Pm  =Pi  +P2  +^3  is  the  total  power 
input  to  the  grid  and  Pe  is  the  total  electric  power  demand. 

2.3.  Electrified  vehicles  and  commuting  patterns 

The  number  of  EVs  is  assumed  to  be  two  million,  which  is  about 
25%  of  the  registered  passenger  vehicles  in  Michigan.  The  commut¬ 
ing  pattern  of  the  vehicle  owners  is  assumed  as  follows:  the  owners 
leave  home  for  work  at  7:20  am  with  2  h  of  standard  deviation;  they 
return  home  around  4:30  pm  with  4.3  h  of  standard  deviation;  the 
battery  state  of  charge  (SOC)  when  plugged-in  is  also  described  by 
a  normal  distribution  with  a  mean  of  0.5  and  a  standard  deviation 
of  0.1 ;  the  vehicles  are  charged  only  during  the  night  (not  at  work). 
The  mean  and  standard  deviation  of  the  commuting  patterns  are 
based  on  the  traffic  count  data  from  Interstate  Highway  5  [20].  The 
battery  capacity  is  assumed  to  be  16  kWh,  all  vehicle  are  assumed 
to  be  charged  by  level-I  chargers  at  1 1 0  V/l  5  A,  and  the  allowed  SOC 
range  is  assumed  to  be  0.3-0.85. 

If  we  simulate  a  scenario  where  every  EV  starts  charging  at  maxi¬ 
mum  power  at  the  moment  they  are  plugged  in,  the  charging  power 
will  add  to  the  grid  load,  as  shown  in  Fig.  2.  The  additional  EV 
demand  increases  the  peak  demand  and  causes  the  peak  load  plants 
to  operate,  which  have  significantly  higher  electricity  generation 
cost. 

2.4.  Renewable  power  source  and  uncertainties 

In  this  study,  electric  power  from  a  wind  farm  is  assumed  to  be 
part  of  the  grid.  The  National  Renewable  Energy  Laboratory  (NREL) 
provides  forecasted  and  observed  wind  power  datasets  for  many 
places  in  the  US  [21].  For  the  algorithm  verification  purposes,  we 
examined  forecasted  and  observed  wind  power  profiles  from  the 
database.  An  example  is  shown  in  Fig.  3.  We  assumed  10%  pene¬ 
tration  rate  of  renewable  power  sources  in  nominal  capacity.  The 
forecasted  profile  will  be  included  in  the  base  load  generation  and 
the  observed  profiles  gives  uncertainties,  the  effect  of  which  is  mit¬ 
igated  by  regulation  services. 


Table  1 

Characteristics  of  the  power  plants  categorized  by  service  types. 


Services 

Base  load 

Load  following 

Regulation 

Control 

Feed  forward 

Integral 

Proportional 

Time  scale 

Hour-day 

10-100  min 

ms-min 

Wind  power  normalized  by  capacity 


Fig.  3.  Example  wind  power  profile  from  the  NREL  database  (June  22nd  2006,  site 
#  3939). 

2.5.  Generation  cost  and  C02  emissions 

A  service  provider  typically  manages  many  power  plants 
that  can  be  categorized  into  base  load  plants,  intermediate  load 
plants,  and  peak  load  plants.  Nuclear  and  coal  power  plants  are 
usually  used  for  the  base  load,  combined  cycle  power  plants  for 
the  intermediate  load,  and  gas  turbines  for  the  peak  load.  In  this 
cost  structure,  we  take  into  account  only  the  energy  cost  so  that 
the  cost  of  regulation  services  is  not  included.  The  generation 
cost  curves  are  based  on  economic  dispatch  rules  for  the  service 
providers.  The  instantaneous  cost  curve  for  the  target  grid  is 
derived  using  the  method  reported  in  [22],  as  shown  in  Fig.  4(a). 
The  carbon  dioxide  emission  curve  is  achieved  by  using  carbon 
dioxide  emission  tables  for  different  fuel  types  [23]  and  plant  types 
[24,25],  as  shown  in  Fig.  4(b).  To  minimize  production  cost  and 
carbon  dioxide  emission,  we  combined  these  two  curves  using  a 
carbon  tax  concept  [26].  The  combined  cost  curve  with  a  carbon 
tax  ($12  per  tonne  of  C02)  is  shown  in  Fig.  5. 

2.6.  Objectives  and  constraints 

The  first  objective  is  to  develop  a  decentralized  algorithm  that 
charges  the  EV’s  battery  as  much  as  possible  while  minimizing  the 
generating  cost  and  carbon  dioxide  emissions  of  power  plants.  The 
secondary  objective  is  to  reduce  the  dependency  of  conventional 
regulation  plants.  To  accomplish  the  objectives,  the  EV  loads  have  to 
be  shifted  to  the  valley  hours  to  minimize  the  operation  of  the  peak 
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Fig.  5.  Combined  cost  curve,  C(P(t)),  (Generation  cost  +  C02  emission),  the  carbon 
tax  of  $12/tC02  is  used. 


load  plants,  and  the  EV  charging  algorithm  should  have  a  frequency 
regulation  capability. 


2.7.  Approaches 

The  overall  grid  system  and  the  charging  controller  are  shown  in 
Fig.  6.  To  achieve  the  design  objectives,  we  developed  a  sub-optimal 
control  algorithm  that  balances  among  multiple  objectives:  min¬ 
imization  of  the  power  generation  cost  and  carbon  dioxide 
emissions.  The  problem  was  first  solved  as  a  linear  programming 
(LP)  problem.  Subsequently,  a  near-optimal  decentralized  control 
algorithm  is  formulated  by  emulating  the  optimal  charging  pattern 
obtained  from  the  LP  solution.  The  optimal  decentralized  control 
algorithm  uses  a  few  pieces  of  information  to  compute  the  charging 
power:  the  forecasted  total  power  demand  of  the  grid,  estimated 
plug  off  time,  and  the  battery  SOC  of  the  local  vehicle. 

For  frequency  regulation,  we  utilize  a  frequency  droop  based 
power  control  algorithm  in  order  not  to  degrade  the  performance 
of  the  load  shifting  algorithm  while  keeping  the  overall  control 
structure  in  a  decentralized  manner. 
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Load  forecasting 
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Charging  Controller 

Load 
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D(l  +  £5) 
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Fig.  6.  The  electricity  grid  model  with  the  proposed  charging  controller. 


3.  Load  shifting 

The  design  process  begins  by  solving  the  LP  optimal  control 
problem  where  all  information  is  available  and  a  central  controller 
determines  the  charging  power  of  each  EV  and  the  power  gener¬ 
ation  of  each  plant.  The  solution  of  the  control  problem  provides 
important  information:  the  optimal  control  pattern  for  each  EV. 
A  decentralized  charging  control  algorithm  can  then  be  obtained 
by  emulating  the  behavior  of  the  LP  solution.  The  remaining  con¬ 
tents  of  this  chapter  include:  formulation  of  the  optimal  control 
problem;  solution  to  the  problem  using  LP  technique;  analysis  of 
the  optimal  control  patterns;  and  derivation  of  the  decentralized 
control  algorithm. 

3.1.  Formulation  of  the  optimal  control  problem 

The  cost  function  to  be  minimized  is  the  electricity  genera¬ 
tion  cost  and  total  carbon  dioxide  emissions.  Given  the  cost  curve, 
C(P(t)),  shown  in  Fig.  5,  the  optimal  control  problem  is  defined  as 
follows: 

T 

min  Vc(P(t)),  (4) 

PEviCn )2-— ' 
t= 0 


where 

N 

P(0  =  rc)  +  PnonEvjt ),  (5) 

n= 1 


C(P(t))  =  Ci  •  min(P(t),  L\ )  +  c2  •  min(max(P(t)  -  Li ,  0),  L2) 

+  c3  •  max(P(t)  -  L\  -  L2,  0),  (6) 

PnonEv(t)  is  the  non-EV  electric  load  and  Pfv(t,  n)  is  the  control  vari¬ 
able  for  t  =  0,  1,  2, . .  .,T,  n  =  1,  2,  3, . .  .,N.  L\  and  L2  are  the  power 
levels  when  the  cost  curve  gets  bended  and  defined  in  Fig.  5.  The 
constraints  for  the  optimization  problem  include 

T 

Y)P£v(t,  n)  ■  At  =  B(n),  (7) 

t= o 

0<PEV(t,n)<PEVUm(t,n),  (8) 


where  B{n)  is  the  total  energy  to  charge  the  battery  of  car  #  n 
as  much  as  it  can  during  the  plug-in  time,  for  n  =  1,  2,  3, . .  .,N,  so 
that  B(n)  is  given  in  this  optimization  formulation.  Pfviim(f)  is  the 
charging  power  limit  of  the  battery,  determined  as  follows: 


f  Pev 

PeV lim(t>  n)  =  < 


when  the  vehicle  n  is  on-line 
when  the  vehicle  n  is  off-line 


(9) 


Eq.  (5)  represents  the  total  power  demand  consisting  of  the  non- 
EV  load  and  the  EV  charging  load;  Eq.  (7)  requires  that  all  EVs  must 
be  charged  as  much  as  unmitigated,  which  prevents  EVs  from  sac¬ 
rificing  charging  completion  for  reducing  power  generation.  Eq.  (8) 
reflects  the  fact  that  the  charging  power  of  each  EV  is  limited  by  the 
maximum  current,  and  no  vehicle-to-grid  power  flow  is  allowed. 
Eq.  (9)  requires  that  EVs  are  charged  only  when  they  are  plugged 


When  the  combined  cost  curve  is  not  linear  then  typically  a 
nonlinear  optimization  problem  with  a  large  number  of  control 
variables  needs  to  be  solved,  which  requires  tremendous  comput¬ 
ing  power.  Fortunately,  the  cost  curve  is  piece-wise  linear  and  has 
increasing  slope.  Therefore,  we  can  use  linear  programming  (LP)  to 
solve  the  problem  by  adding  more  control  variables  and  constraints. 
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3.2.  Solution  through  linear  programming 

The  formulated  optimization  problem  can  be  solved  by  modify¬ 
ing  the  cost  function,  as  follows: 

T 

min  (ci  •  9i(0  +  c2  •  92(f)  +  c3  •  q3(t)),  (10) 

where  0<Ci  <c2<c3  are  defined  in  Fig.  5.  The  constraints  consist 
of  linear  constraints  and  conditional  constraints.  The  linear  con¬ 
straints  are 


N 


9l(f)  +  92(0  +  93(0  =  y  PEvit,  n)  +  PnonEv{t), 
n= 1 

(id 

T 

yPEv(t ,  n)  =  B(n), 

(12) 

t= 0 

0<P£V(t,n)<P£Vlim(t,n), 

(13) 

0  <  9i(f)  <  Li ,  0  <  q2(f)  <  L2,  0  <93(t). 

(14) 

And  the  conditional  constraints  are 


(a)  Optimal  charging  power 


1  2 
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Time  /  hour 


(b)  Power  demand 


(c)  Example  charging 
power  profiles 


j*: 


Fig.  7.  Optimal  results  obtained  from  LP:  (a)  the  charging  pattern  of  EVs;  (b)  the 
resulting  electric  power  demand;  and  (c)  three  examples  of  charging  patterns  for 
different  SOC  levels. 


. _  .  r  f  <ji(t)  =  p(t), 

if  P(t)<Lu  < 

l  <?2(t)  =  0,  <J3(t)  =  o, 


if L,  <  P(t)  <  L2, 
ifP(t)  >  Lj  +L2, 


1 


f  =  Lj, 

q2(t)  =  P(t)-L1,  q3(t)  =  0, 
f  qj(t)  =  Lu  q2(t)  =  L2 , 

l  <J3(t)  =  P(t)-Li -L2, 


3.4.  Development  of  the  decentralized  control  algorithm 

To  analyze  the  charging  pattern,  the  EVs  are  grouped  by  the 
plug-off  time,  as  shown  in  Fig.  8.  The  charging  power  of  EVs  that 
(15)  have  the  same  plug-off  time  shows  linear  dependency  to  SOC,  and 
the  sensitivity  changes  over  time.  As  a  result,  the  charging  power 
can  be  expressed  as: 


PEv,SFr(t ,  n)  =  I<(t,  n)  •  (SOCmax  -  SOC(t ,  n)), 


(16) 


where  qi(t),  q2(t),  and  q3(t)  vary  only  in  the  linear  ranges  of  the 
piece-wise  linear  cost  curve.  The  introduction  of  new  variables 
transforms  the  nonlinear  cost  function  to  a  linear  cost  function  and 
in  the  meantime  adds  the  conditional  constraints.  Flowever,  these 
conditional  constraints  are  not  active  (the  optimal  solution  does  not 
exist  in  the  hyper-plane  defined  by  the  conditional  constraints), 
as  proved  in  Appendix.  Therefore,  we  can  ignore  the  conditional 
constraints  and  the  problem  can  be  solved  using  a  standard  linear 
programming  technique. 


where  SOCmax  is  the  maximum  allowable  SOC  that  is  0.85  in  this 
problem.  The  gradient  JC(t,  n)  is  computed  from  the  LP  solutions 
and  plotted  in  Fig.  9.  The  gradient  increases  as  the  time  to  plug-off 
decreases,  and  as  the  time  to  valley-end  decreases.  Also,  the  amount 
of  available  cheap  power  affects  I<(t,  n). 

Another  way  to  derive  Pr/seKC  n)  is  based  on  the  observation 
that  charging  power  is  highly  related  to  the  amount  of  available 
cheap  power,  which  is  defined  as  L\  +  L2  -  PnonEv{t).  For  perfect 
valley-filling,  the  following  equation  must  hold: 


3.3.  Optimal  results  and  analysis 

The  parameters  of  the  optimization  problem  are  as  follows: 
T=24  with  1  h  step,  N=  100  (the  two  million  EVs  are  grouped  into 
100  sets).  The  plug-in  time,  plug-off  time,  and  the  initial  SOCs  are 
randomly  generated  based  on  the  commute  pattern,  and  their  val¬ 
ues  are  assumed  to  be  known  to  the  centralized  controller.  The 
problem  is  solved  using  the  LP  solver  in  MATLAB  (linprog)  and  the 
results  are  plotted  in  Fig.  7.  In  the  optimal  solution,  the  additional 
load  due  to  EV  charging  fills  the  valley  of  the  demand  profile  during 
the  night  so  that  the  peak  power  plant  does  not  need  to  operate.  If 
we  compare  the  demand  profile  of  Fig.  7(b)  with  that  of  Fig.  2,  we 
can  see  that  the  EV  charging  demand  is  shifted  to  early  morning 
hours.  Fig.  7(c)  shows  examples  of  optimal  charging  patterns.  The 
EVs  begin  to  charge  as  soon  as  the  valley  period  begins,  which  is 
defined  as  the  time  when  the  base  demand  is  below  the  peak  load 
line,  where  we  do  not  need  to  turn  on  the  peak  load  plants.  The  EVs 
also  finish  their  charging  when  the  valley  period  ends.  The  charging 
power  depends  on  the  battery  SOC:  the  EV  with  a  lower  SOC  uses 
higher  charging  power.  The  charging  power  also  depends  on  the 
amount  of  available  grid  power. 
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^PEVMU  n)  =  l\  +  L2  -  PnonEv(t)- 
n= 1 


(17) 
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Fig.  8.  Optimal  charging  power  from  the  LP  solutions  was  found  to  be  proportional 
to  “SOC  deficit”. 
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(a)  3D  plot  (b)  Contour  plot 


Then,  the  average  charging  power  per  vehicle  at  time  t  is  calculated 
as  follows: 


n  rn  _  Yln=\pEV,SFr{t,  n)  _  Li  +L2  —  PnonEv{t) 

EV'SFT[  )  ~  N^Ki)  “ 


(18) 


where  NEv{t)  is  the  number  of  EVs  that  are  plugged  in.  In  this  stage, 
PnonEvit)  and  NevW  are  assumed  known  for  all  t.  Using  this  average 
charging  power,  we  propose  another  charging  algorithm: 


PEV,SFr(t ,  n)  =  R(t,  n)  •  PEv,SFr(t)> 


(19) 


where  R(t,n)  is  the  charging  gain.  Because  the  EV’s  battery  should 
be  fully  charged  by  the  end  of  valley  hours  or  the  plug-off  time,  the 
following  equation  must  hold: 

pTn 


socn 


-soc(,'n)-cfe/ 


Pfv(r,  n)  •  dr 


—  [Tn  i 


=  /  Rfa  n)  •  Pev,sftM  •  dr, 


(20) 


where  Tn  =  min{tplug.offtn,  tvaUey_end)  and  C(n)  is  the  battery  capacity 
of  EV  #n.  Because  we  do  not  know  future  R( r,  n),  it  is  assumed  to 
be  time  invariant  and  we  have: 


=  C(n)(SOCmax-SOC(t,n)) 

flend  Pev,sft{t)  •  dr 


This  calculated  R(t,  n)  is  assumed  to  be  constant  until  PEv,SFr(t) 
and  SOC{tf  n )  are  updated.  Plug  (21)  into  (19),  we  have 


PEv,SFr(t,  n)  =  ^n_)  'P£^sfT(t)  (SQCmax  _  S0C(t,  n));  (22) 

ft  n  Pev,setM  ■  dr 

which  is  the  equation  of  the  decentralized  charging  algorithm.  We 
should  point  out  that  K(t,  n )  in  Eq.  (16)  can  be  calculated  by  com¬ 
paring  Eqs.  (16)  and  (22): 


K(t,n)=  y_)  -*W0  , 
ft  n  Pev,sftM  •  dr 


(23) 


where  Pev,SFT (0  =  P\  +  ^2  -  PnonEv{t)/^Ev{t)- 

The  gradient  JC(t,  n)  in  Eq.  (23)  should  be  verified  by  comparing 
against  the  LP  solution.  To  see  the  robustness  of  this  gradient  func¬ 
tion,  two  different  power  profiles  are  used.  Fig.  10  compares  the 
gradient  values  from  the  optimization  results  and  from  Eq.  (23). 
An  imagined  double  valley  profile  is  used  as  a  challenging  case. 
Even  though  the  gradient  function  was  derived  using  a  single  valley 
profile,  it  matches  well  with  LP  results  calculated  from  the  double 
valley  profile,  which  shows  the  robustness  of  the  gradient  func¬ 
tion.  The  parameters  of  the  controller  should  be  determined  by  two 
separate  players.  C(n)  is  determined  by  the  EV  designer  and  L\  L2 
and  NEv{t)  should  be  by  the  grid  operator  in  consideration  of  the 
generation  cost  structure. 


5 


o 


"O 

c 

03 

E 

0 

0 


(a)  Single  valley  case 


LP 

Function 


Time  /  hour 


Time  /  hour 


(b)  Double  valley  case 


Time  /  hour 


Fig.  10.  Verification  of  the  robustness  of  the  gradient  function  K(t,  n). 


The  control  algorithm  in  Eq.  (22)  can  be  implemented  in  a 
decentralized  fashion,  which  consists  of  two  parts:  gain  and  SOC 
deficiency.  The  gain  becomes  larger  when  the  time  to  plug-off 
is  shorter,  and  when  the  amount  of  cheap  power  is  higher.  The 
dependency  on  SOC  deficiency  ensures  that  vehicles  with  lower 
SOC  receive  higher  charging  power.  Additionally,  the  dependencies 
show  simple  proportional  relations.  To  implement  this  algorithm, 
predicted  future  demand  from  current  time  to  the  end  of  the  valley 
hours,  and  the  plug-off  time  of  the  vehicle  are  required.  An  accu¬ 
rate  forecast  of  non-EV  electric  load  is  already  available  and  used 
in  the  power  transmission  industry,  therefore,  the  end  of  the  prob¬ 
lem  horizon  (end  of  valley  hours)  is  readily  available.  The  plug-off 
time  of  individual  EV  needs  to  be  known,  perhaps  estimated  from 
user  input  or  learned  from  past  behaviors.  Also,  we  assume  that  the 
number  of  EVs  connected  to  the  grid  can  be  estimated  from  histor¬ 
ical  data,  or  from  the  binary  on-off  data  from  their  smart  chargers. 
Therefore,  the  controller  is  implementable  in  a  decentralized  way 
using  today’s  technology.  The  control  algorithm  is  simple  and  the 
computation  can  be  done  locally.  Information  from  vehicle  to  grid 
is  not  required— except  the  binary  plug-in  data.  The  only  informa¬ 
tion  that  needs  to  broadcast  from  the  grid  to  the  vehicles  is  the 
predicted  cheap  power  trajectory. 

4.  Frequency  regulation 

Another  important  attribute  of  EV  charging  is  its  fast  response 
time.  EV  batteries  have  a  faster  power  slew  rate  than  conventional 
power  generators  such  as  steam  turbines.  Due  to  this  characteris¬ 
tic,  EVs  can  replace  some  of  the  regulation  service  units.  Battery 
charging  will  not  be  deteriorated  by  the  regulation  control  because 
regulation  is  a  zero-energy  service  that  compensates  for  minute- 
to-minute  fluctuations  [3].  Long-horizon  services  (>1 0  min)  such  as 
load  following  or  energy  balancing  will  continue  to  be  provided  by 
other  power  sources. 

The  regulation  control  using  EV  batteries  is  designed  as  a  fre¬ 
quency  droop  control: 


Pev,reg( t,  n)  =  n )  -  co0),  (24) 

where  co(t,  n)  is  the  frequency  measured  by  EV  #  n  at  time  t.  kj{t ) 
is  the  droop  control  gain  that  is  time  varying  and  is  broadcast  by 
the  grid  operator.  From  the  point  of  view  of  the  grid,  the  effective 
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Fig.  11.  The  final  control  structure.  Grid  operator  broadcast  signals  to  each  EV.  The  controller  inside  the  dashed  box  is  located  in  each  EV. 


gain  ( keff )  depends  on  the  number  of  EVs  connected  to  the  grid,  as 
follows: 

N 

PregW  =  'y^PEV,REG(t,  n)  =  NEv{t)kf{t){o){t ,  n)  -  co0) 

n= 1 

=  n)  -  coq).  (25) 

A  constant  effective  gain  is  preferable  for  consistent  regulation 
performance.  However,  when  NEV(t )  is  very  small,  constant  effec¬ 
tive  gain  cannot  be  achieved.  Therefore,  we  propose  the  following 
adaptive  gain: 

k'(t)=N^r  (26) 

where  Kf  is  the  desired  overall  gain  and  NEViNOmW  is  the  nominal 
number  of  vehicles  that  only  need  to  be  updated  periodically.  I<f  is 
complementary  to  kp  in  Eq.  (2)  and  they  determine  the  amount  of 
the  frequency  droop.  The  control  parameters  of  regulation  plants, 
such  as  I<f ,  kp,  are  usually  determined  from  the  maximum  available 
regulation  power  dividend  by  the  maximum  allowable  frequency 
droop.  Mismatch  between  the  true  number  of  EVs  and  its  nominal 
value  will  cause  slight  variation  in  effective  regulation  gain.  When 
the  total  number  of  EVs  is  small,  they  may  not  provide  adequate 
regulation  capability. 

5.  Algorithm  integration 

The  load  shifting  algorithm  and  the  regulation  algorithm  were 
designed  to  achieve  their  individual  objectives.  However,  in  con¬ 
cept,  the  algorithms  will  not  interfere  with  each  other  because 
of  their  time-scale  separation.  Therefore,  the  two  controllers  are 
combined  in  a  simple  way: 

PEv(t,  n)  =  sat{PEV,sFr(t,  n)  +  PEv,REc(t,  n), 

~  PeV lim(0>  PeV lim(0>  }>  (27) 

where  sat  is  a  saturation  function  that  takes  into  account  the  power 
limit  of  EV  batteries.  The  overall  control  structure  is  shown  in 
Fig.  11.  The  controller  in  the  grid  operator  side  computes  PEv,sft (0 


and  kj(t)  to  broadcast  to  EVs.  The  broadcasting  is  done  at  a  very  low 
time  frequency  such  as  every  ten  minutes.  The  controller  in  the  EV 
charger  computes  charging  power  for  load  shifting  and  for  regula¬ 
tion  using  local  measurement  signals:  SOC  and  AC  frequency.  This 
controller  performs  computation  at  a  higher  rate  such  as  every  10  s. 

6.  Algorithm  verification 

The  decentralized  control  algorithm  was  verified  using  six 
simulation  scenarios  categorized  into  three  groups,  as  summarized 
in  Table  2.  The  first  group  was  designed  to  verify  the  performance 
of  the  load  shifting  algorithm  under  perfect  information.  Scenario 
1  uses  the  demand  profile  shown  in  Fig.  1.  Scenario  2  uses  a 
demand  profile  that  has  double  valleys  shown  in  Fig.  10(b).  The 
second  group  was  designed  to  verify  the  performance  of  the  load 
shifting  controller  under  imperfect  estimation.  Scenario  3  uses  the 
single  valley  profile  with  inaccurate  NEV(t):  the  estimated  NEV(t)  is 
assumed  to  be  delayed  by  1  h  from  the  true  values.  Scenario  4  also 
uses  the  same  demand  profile  but  the  total  number  of  EVs  is  50% 
higher  than  the  expected  (2  million  vehicles  were  expected  but  3 
million  vehicles  were  plugged  in).  The  last  group  was  designed  to 
evaluate  the  integrated  algorithm  (load  shifting  and  regulation). 
Scenario  5  uses  the  single  valley  demand  profile  with  a  wind 
power  source.  Scenario  6  uses  the  same  demand  profile  and  the 
wind  power  source  but  the  regulation  controller  was  turned  off 
to  investigate  the  effect  of  EV’s  regulation  operation.  In  these 
scenarios,  we  assumed  perfect  estimation  except  wind  power 
forecasting.  The  rest  of  the  simulation  conditions  were  common 
for  every  scenario  as  follows:  the  simulation  horizon  was  24 h; 
the  time  step  was  10  s;  the  plug-in/off  time  and  the  initial  SOCs 
were  randomly  generated  but  identical  for  each  simulation;  and, 
the  number  of  vehicle  fleets  was  1000,  which  means  that  each  one 
represented  2000  EVs  (2500  EVs  for  Scenario  4). 

6. 1 .  Load  shifting  algorithm  with  perfect  estimation 

The  simulation  results  for  Scenarios  1  and  2  are  plotted  in 
Fig.  12  and  the  performance  summary  is  shown  in  Table  3.  In  both 
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Table  2 

Simulation  scenarios  to  verify  the  performance  of  the  proposed  control  algorithms. 


Scenario 

Controller 

NEv(t)  estimation 

Renewable  power 

Demand  profile 

1 

Load  shifting 

Perfect 

No 

Single  valley 

2 

Load  shifting 

Perfect 

No 

Double  valley 

3 

Load  shifting 

1  h  delayed 

No 

Single  valley 

4 

Load  shifting 

25%  less 

No 

Single  valley 

5 

Shifting  +  regulation 

Perfect 

Yes 

Single  valley 

6 

Load  shifting 

Perfect 

Yes 

Single  valley 

scenarios,  the  baseline  controller  allows  unmitigated  charging. 
Therefore,  it  achieved  the  best  battery  charging  performance  but 
the  power  generation  cost  was  much  higher  than  the  controlled 
cases.  The  decentralized  load  shifting  algorithm  effectively  shifted 
the  peak  load  into  the  valley,  achieving  lower  generation  cost  (5-8% 
reduction)  without  losing  the  battery  charging  performance. 


6.2.  Load  shifting  algorithm  with  imperfect  estimation 

Scenario  3  was  designed  to  test  the  robustness  against  num¬ 
ber  of  electrified  vehicles,  NEV(t).  We  assumed  that  the  estimated 
vehicle  number  was  a  time  delayed  version  of  the  actual  num¬ 
ber  and  the  results  are  shown  in  Fig.  13(a)  and  (b).  Even  though 
the  inaccurate  estimation  slightly  degraded  the  charging  perfor¬ 
mance  (from  99.99%  of  Scenario  1  to  99.93%  of  Scenario  3)  the 
valley  filling  worked  well  and  the  degradation  was  small.  Scenario 
4  was  designed  to  study  the  effect  of  more  EVs  than  expected  and 
the  results  are  shown  in  Fig.  13(c)  and  (d).  The  expected  num¬ 
ber  of  EVs  was  2  million  but  the  number  of  connected  EVs  was 
3  million.  The  required  charging  energy  was  much  higher  than  the 
energy  available  in  the  valley  hours.  As  a  result,  the  more  expen¬ 
sive  power  source  is  used.  However,  battery  charging  still  happens 
in  an  orderly  manner  and  the  charging  performance  was  still  high 
(99.92%),  as  shown  in  Table  3.  These  two  simulations  show  that  the 
performance  of  the  decentralized  load  shifting  algorithm  is  quite 
robust  and  vehicle  charging  completion  took  priority  over  other 
objectives.  To  achieve  leveled  valley  filling  in  the  case  of  3  million 
EVs,  the  grid  operator  has  to  compute  Pev,sft(0  using  larger  L\  +L2 
than  in  the  case  of  2  million  EVs. 


6.3.  Integrated  controller  with  wind  power  uncertainties 

Simulations  with  the  integrated  controller  were  performed  to 
evaluate  two  important  aspects:  whether  the  load  shifting  per¬ 
formance  and  the  final  SOC  are  degraded  by  integration  of  the 
regulation  controller  and  how  the  regulation  controller  is  ben¬ 
eficial.  Scenario  5  uses  the  integrated  (load  shifting  +  regulation) 
algorithm  with  the  single  valley  demand  profile  and  the  wind 
power  source  profile  shown  in  Fig.  3.  Scenario  6  uses  only  the 
load  shifting  controller  with  the  single  valley  demand  profile 
and  the  wind  power  source.  In  both  scenarios,  the  forecasted 
demand  profile  and  forecasted  wind  power  profile  were  used  to 
determine  the  required  power  generation  from  the  base  plants. 
The  wind  power  uncertainties  are  compensated  by  the  regula¬ 
tion  services  from  the  conventional  regulation  plants  and  EVs 
in  Scenario  5  and  only  by  the  conventional  regulation  plants  in 
Scenario  6. 

The  results  are  plotted  in  Fig.  14  and  summarized  in  Table  4. 
The  load  shifting  performance  and  the  charging  performance  of 
Scenario  5  were  not  degraded  compared  to  those  in  Scenario  6, 
which  confirms  that  the  regulation  controller  does  not  degrade  the 
load  shifting  performance  and  the  charging  performance,  which 
is  expected.  Furthermore,  the  regulation  controller  improved  the 
frequency  regulation  with  reduced  reliance  on  conventional  reg¬ 
ulation  power  plants,  which  indicates  the  possibility  of  cost 
reduction  in  regulation  services. 

In  summary,  as  shown  in  Fig.  15,  the  decentralized  charging 
control  algorithm  reduces  generation  cost  and  C02  emissions  with¬ 
out  degrading  EV  charging  completion  rates.  In  the  meantime,  the 
algorithm  contributes  to  the  frequency  regulation  services.  The 
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Fig.  12.  Verification  of  the  load  shifting  controller  under  perfect  estimation. 
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Fig.  13.  Verification  of  the  load  shifting  controller  under  imperfect:  (a)  when  NEV(t ) 
estimation  is  inaccurate,  (b)  number  of  total  EVs  are  more  than  expected. 
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Table  3 

Performance  of  the  load  shifting  algorithm  compared  to  the  baseline  and  LP  optimal  results. 
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Scenario 

Controller 

Generation  costa 

CO2  emission3 

Charging  performance13  (%) 

Baseline 

1 

1 

100 

1 

LP  optimal 

0.907 

0.969 

100 

Decentralized 

0.923 

0.976 

99.99 

Baseline 

1 

1 

100 

2 

LP  optimal 

0.908 

0.970 

100 

Decentralized 

0.948 

0.985 

99.98 

3 

Baseline 

1 

1 

100 

Decentralized 

0.923 

0.976 

99.93 

4 

Baseline 

1 

1 

100 

Decentralized 

0.925 

0.976 

99.92 

a  Generation  cost  and  C02  emission  are  normalized  by  the  values  of  the  baseline  cases. 

b  Charging  performance  is  the  sum  of  all  served  SOCs  divided  by  the  sum  of  all  SOCs  of  the  baseline  controller  case. 


Table  4 

Performance  of  the  integrated  algorithm  compared  to  the  load  shifting  algorithm  only  case. 


Scenario 

Controller 

Generation  cost 

Charging 
performance  (%) 

Regulation 
power3  (MW) 

Frequency 
deviation13  (Hz) 

5 

Shifting  +  Regulation 

0.918 

99.98 

0.116 

0.069 

6 

Load  shifting 

0.917 

99.98 

0.146 

0.087 

a  Regulation  power  is  the  root-mean-squared  power  output  of  the  conventional  regulation  plant. 
b  Frequency  deviation  is  the  root-mean-squared  frequency  deviation  from  60  Hz. 
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Fig.  14.  Verification  of  the  integrated  algorithm:  the  performance  of  the  integrated 
algorithm  compared  to  the  load  shifting  algorithm  only  case. 


load  shifting  performances  of  decentralized  controller  are  closed 
to  the  LP  optimal  performances  and  the  frequency  regulation  per¬ 
formance  is  four  times  better  than  the  unmitigated  case  while 
reducing  regulation  power  by  75%. 


7.  Conclusions 

This  paper  presents  the  design  of  a  decentralized  charging 
algorithm  for  electrified  vehicles.  The  algorithm  had  two  objec¬ 
tives:  load  shifting  for  efficiency  and  frequency  regulation  for 
power  quality.  The  load  shifting  algorithm  mimics  the  behavior 
of  a  global  optimal  solution  obtained  through  the  linear  program¬ 
ming  technique.  The  developed  load  shifting  algorithm  requires 
four  pieces  of  information:  forecasted  base  load  profile,  estimated 
number  of  plugged  vehicles,  estimated  plug-off  time,  and  the 
battery  SOC  of  the  vehicle  being  charged.  The  regulation  algo¬ 
rithm  is  based  on  power-frequency  droop.  Simulation  studies  show 
that  the  proposed  algorithm  minimizes  electricity  generation  cost 
and  C02  emissions  and  reduces  the  usage  of  the  conventional 
regulation  power  plants  without  compromising  battery  charging 
performance. 


Cost  Emission  Performance  Power  Deviation 


Fig.  15.  Comparison  of  performance  indices  normalized  by  the  unmitigated  case 
(Scenario  5).  The  performance  of  frequency  regulation  is  not  applicable  to  the  LP 
optimal  results. 
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Appendix. 

Suppose  an  optimal  solution  exists  when  the  conditional  con¬ 
straints  shown  in  Eq.  ( 1 5 )  are  violated  at  t\ .  Note  that  the  following 
conditions  hold  at  any  time. 

P(t i )  =  <?i(ti )  +  q2(t i )  +  q3(ti ),  0  <  <Ji(t)  <  U , 

0  <  Q2(t)  <  L2,  0  <  q3(t),  0<c1<c2<c3. 
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(i)0  <P(ti)<L, 


ti-i 


J*  C(P(0)  +  Ci <?,(t1)  +  C2Cj2(ti)  +  c3q3(ti) 

t=0  t=t!+l 

t"!  —  1  T 

=  £c(P(t))+  ^  C(P(t))  +  c1P(t1)-c1P(t1)  +  c1q1(t1)  +  c2q2(t1)  +  c3q3(t1) 

t=0  t=t]+l 

t,-l  T 

=  y^C(P(t))  +  ^  C(P(t))  +  c,P(t! )  -  c, (q,(t! )  +  q2(ti )  +  Q3(ti ))  +  Cjq, (t, )  +  c2q2(h )  +  C3<73(ti ) 


t= 0  t=t!+l 

tl-1  T  q-1  T 

=  ^C(P(t))  +  c(p(f))  +  CiP(ti)  +  (c2-c1)q2(t1)  +  (c3-c1)q3(ti)>^C(P(t))+  ^  C(P(t))  +  dP^ ) 

t=0  t=t!+l  t=0  t=ti+ 1 


The  equality  holds  only  when  ^2(^1 )  =  <^3(^1 )  =  0*  There¬ 
fore,  the  optimal  solution  is  achieved  when  qi(ti)  =  P(ti)  and 
<72(^1 )  =  <73(ti)  =  0- 
(ii)  U  <P(fi)<^i  +L2 


J* 


q- 1  T 

=  £c(P(t))+  ^  C(P(f))  +  c1q1(t1)  +  c2q2(t1)  +  c3q3(t1) 

t=0  t=q+l 

q-1  T 

=  yC(P(t»  +  C(P(t))  +  CiL!  +C2L2  +c3(P(ti)-Li  -L2)-CiL!  -  C2L2  -c3(P(ti)-L!  -  T2)  +  Ci<7i(£i )  +  C2<?2(fi )  +  C3q3(ti ) 

t=0  t=q  + 1 

q-l  T 

=  yC(P(t))  +  C(P(t))  +  CiLi  +C2L2  +c3(P(ti)-Li  -P2)  +  (c3  -Ci)(Li  -qi(ti))  +  (c3  -  c2)(L2  -^2(^1)) 

t=0  t=q+l 

q-1  T 

>y^C(P(t))+  ^  C(P(t))  +  c1L1+c2L2  +  c3(P(ti)-L1-L2) 
t=o  t=q+i 


The  equality  holds  only  when  qi(ti)  =  Li,  q3(ti)  =  0.  There¬ 
fore,  the  optimal  solution  is  achieved  when  <h(ti)  =  Ii, 
q2(t1)  =  P(ti)-L1,and  q3(ti)  =  0. 

(iii)P(t1)>L1+L2 


J*  =  yW))+  ^^C(P(t))  +  c,qi(t1)  +  c2q2(ti)  +  C3q3(ti) 

t=0  t=t1  +1 

tl-1  T 

=  y^C(P(t))+  yy  C(P(t))  +  CiLi  +C2I2  +c3(P(ti)-Ii  -I2)-Cili  -C2I2  -c3(P(ti)-Ii  -I2)  +  Ciqi(ti)  +  C2q2(ti)+C3q3(ti) 

t=0  t=ti+l 

t!-i  r  ti  —i  r 

=  ^C(P(t))+^C(P(t))+c1L1+c2L2+C3(P(t1)-L1-L2)+(C3-c1)(L1-q1(t1))+(C3-C2)(L2-q2(t,))>^C(P(t))+  ^C(P(t))  +  c,L,  +  c2L2  +  c3(P(t,)  -  L,  -  L2) 


The  equality  holds  only  when  <ji(ti)  =  Ii,  q2{h)  =  L2.  Therefore, 
the  optimal  solution  is  achieved  when  qi(ti)  =  Li,  q2ih)  =  L2l  and 

^3  (  f  1 )  =  P(  f  1 )  —  ^1  —  f-2* 

The  results  of  (i)-(iii)  contradict  the  assumption,  which  implies 
that  the  conditional  constraints  are  satisfied  when  the  optimal  solu¬ 
tion  exists,  in  other  words,  the  constraints  are  not  active. 
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